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ABSTRACT 

It is widely recognized that cold dark matter models predict abundant dark 
matter substructure in halos of all sizes. Galaxy-galaxy lensing provides a unique 
opportunity to directly measure the presence and the mass of such substructures 
in clusters of galaxies. Here we present the mass function of substructures ob- 
tained from lensing in massive HST cluster lenses, and compare it with that 
obtained from high-resolution cosmological N-body simulations. We find excel- 
lent agreement in the slope and amplitude in the mass range 10 11 — 10 12 ' 5 M 
probed by the observations, highlighting a significant success of the CDM theory. 
At lower mass scales, the simulations predict a large abundance of small sub- 
structures below our detection threshold with galaxy-galaxy lensing, with many 
of them being plausibly associated with faint cluster galaxies. Our results sug- 
gest that the CDM substructure abundance on the scale of clusters is in good 
agreement with present observational data. 

Subject headings: cosmology: dark matter — gravitation: gravitational lensing 

1. Introduction 

Arguably the strongest challenges presently faced by the otherwise highly successful cold 
dark matter (CDM) paradigm for structure formation in the Universe occur on 'small scales', 
i.e. within individual dark matter halos. One of the most prominent of these problems is the 
long standing and controversially debated question whether rotation curves of low surface 
brightness galaxies can be accommodated in the CDM theory (e.g. McGaugh & de Blok 
1998; Hayashi et al. 2003; Navarro et al. 2004). Another important issue is concerned with 
the amount of dark matter substructure theoretically predicted by CDM, and whether this 
is in conflict with the relative paucity of satellites observed in the Milky Way and the Local 
Group (e.g. Moore et al. 1999; Klypin et al. 1999; Stoehr et al. 2002). 
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Interestingly, both of these problems have only been recognized once N-body simula- 
tions reached high enough resolution for detailed studies of the phase-space structure of 
dark matter halos, allowing them to probe the inner dark matter cusp to sub-kpc scales, 
and to establish that dark halos are filled with a host of self-bound dark matter satellites, 
instead of being the smooth halos envisioned based on earlier low-resolution work. It has 
now become clear that the existence of substructure is a generic prediction of hierarchical 
structure formation in CDM models, where the assembly of collapsed mass proceeds via a 
merger hierarchy that allows some of the infalling dense lumps to survive as dynamically 
distinct substructure inside virialized halos until late times. 

On the scales of ordinary galaxies like the Milky Way (~ 1O 12 M ), the number of 
bound satellites predicted by the simulations is much larger than the dozen or so satellites 
actually detected around our Galaxy. Comparing the observed velocity distribution function 
of these satellites with the predicted one (Moore et al. 1999; Klypin et al. 1999), this has been 
interpreted as evidence for a serious surfeit of dark matter satellites, and even considered 
to mark a "crisis" of the ACDM model. However, as Stoehr et al. (2002) have shown, 
there is significant uncertainty in the conversion of stellar line-of-sight velocity dispersion 
into peak circular velocities of dark matter substructures. Using high-resolution simulations 
of a Milky Way sized halo, Stoehr et al. in fact were able to show that all the known 
satellites of the Milky Way can be comfortably hosted by the most massive dark matter 
substructures expected for its halo. In this picture, there is however still a sea of small 
mass dark matter sub-halos which are largely devoid of stars. The leading hypothesis to 
explain this situation is to allude to baryonic processes of galaxy formation, inhibiting star 
formation preferentially in low mass halos. Proposed processes for such feedback include a 
photo-ionizing UV background or the expulsion of gas in shallow potential wells by supernova 
explosions (e.g. Bullock et al. 2000; Benson et al. 2002; Kravtsov et al. 2004). 

If vast numbers of dark sub-halos exist, lensing might be the best way to detect them 
(e.g. Trentham et al. 2001). This is not only true for the scales of galaxies, but particularly 
for rich clusters of galaxies, where strong and weak lensing effects can be combined in a 
powerful way to construct detailed mass models for clusters. This allows in principle a direct 
mapping of the sub-halo mass function, thereby providing an important test of the CDM 
paradigm. Upon comparison with the observed galaxy luminosity function of clusters, it can 
also give valuable insights into the galaxy formation process. 

In this letter, we present results of the first comparison between the lensing determined 
substructure mass function in clusters on mass scales ranging from about 10 11 — 10 12 ' 5 M 
with that inferred from high resolution cosmological N-body simulations. To this end, we 
construct high resolution mass maps of galaxy clusters by applying galaxy-galaxy lensing 
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techniques where substructure inside clusters is detected and mapped using the anisotropics 
that they produce in the observed shear field. We compare our lensing measurements directly 
with results from high-resolution N-body simulations, allowing us to test the robustness of 
the CDM model and the associated hierarchical galaxy formation paradigm. Note that 
compared to the scales of galaxies, we expect many more dark matter structures to be 
visible optically in clusters, making the comparison of sub-halo mass functions less affected 
by uncertainties in the galaxy formation physics. Therefore, full consistency between the 
abundance of optically detected galaxies, substructure in CDM N-body models, and sub- 
halos detected by lensing can be asked for; establishing such a consistency can be viewed as 
a strong test of the theoretical paradigm. Below, we outline our methodology for mapping 
substructure in clusters. We then present the mass function derived from the maximum- 
likelihood analysis. This is followed by a comparison with the substructure detected in 
simulated clusters, and a discussion of our results. We adopt h = 0.7, f2 — 0.3 and 
tt A = 0.7. 



2. Quantifying substructure in clusters with gravitational lensing 

We obtain the mass spectrum of clumps in a cluster by combining constraints from 
strong lensing observations (highly magnified multiply imaged systems) and weak lensing 
(the radially averaged tangential shear profile). To this end, we use a self-similar, parametric 
mass model to describe the cluster as a composite of a large-scale smooth mass distribution 
and several sub-clumps (Natarajan & Kneib 1997; Natarajan et al. 1998, 2002, 2004). These 
sub-halos are associated with bright, early-type galaxies in the cluster under the assumption 
that mass traces light. The local anisotropies in the shear field induced by the sub-halos 
in their vicinity is then used statistically to quantify the mass of the sub-clumps and their 
spatial extents. A likelihood method is used to retrieve characteristic halo parameters (e.g. 
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Natarajan & Kneib 1997; Natarajan et al. 1998; Geiger & Schneider 1998). On applying 
these techniques to an ensemble of HST cluster lenses (results are presented in Table 1) we 
find that the spatial extents inferred are consistent with tidal stripping; early-type galaxies 
do possess dark halos that extend well beyond the light but these halos are more compact 
than those around field galaxies of equivalent luminosity. 

In performing the likelihood analysis to obtain characteristic parameters for the sub- 
clumps in the cluster we assume that light traces mass. This is an assumption that is well 
supported by galaxy-galaxy lensing studies in the field (Wilson et al. 2001) as well as in 
clusters (Clowe & Schneider 2002; Hoekstra 2003). The individual galaxies and the smooth 
cluster component are modeled self-similarly with truncated pseudo-isothermal mass distri- 
butions (PIEMD). The parameters that characterize a truncated PIEMD are: a truncation 
radius r t identified with the tidal radius, a core radius r and a central velocity dispersion 
o"o. For the smooth component the values of these parameters are set by the observed strong 
lensing features, and for the galaxies, combined constraints from the strong lensing and the 
weak shear field determine the best-fit parameters for a fiducial galaxy halo. These values 
recovered from the maximum likelihood analysis are shown in Table 1. In order to relate 
the observed light distribution in the early-type cluster galaxies to their masses a set of 
physically motivated scaling laws are assumed: 

a ° = ^ ; r ° = r °* (£) ; n = rt * (zi) • (1) 

The total mass of the sub-halo associated with a galaxy of luminosity L is: 

„ 2 ( L\ 1/2+a M 2 fL\ 1/2 - a 

M(xa °* rt *{r*) ; t K a °* Tt * [ii ) ■ (2) 

Note however that for our choice of mass model and value adopted for the exponent a, the 
mass to light ratio is not constant with radius within an individual galaxy halo. The derived 
mass spectrum of sub-halos is not a strong function of a, as discussed in (Natarajan et al. 
2004). 

Dark halos are associated with the locations of bright, early-type cluster galaxies and 
the fiducial parameters for a typical halo are then extracted from the likelihood analysis. 
A high resolution mass model for the entire cluster is built using the strong lensing regime 
to constrain the inner region (r ~ ^Einstein), and the local anisotropies in the shear field 
are used to obtain properties of the galaxy halos around early-type cluster members. Since 
the procedure involves a scaled, self-similar mass model that is parametric, we obtain a 
mass estimate for the dark halos (sub-clumps) of the cluster galaxies as a function of their 
luminosity. This provides us with a clump mass spectrum. Note that tidal truncation 
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by the cluster causes these halo masses to be lower than that of equivalent field galaxies 
at comparable redshifts obtained from galaxy-galaxy lensing. The fraction of mass in the 
clumps is only 10-20% of the total mass of the cluster within the inner 500 h^ 1 kpc of these 
high central density clusters. We are limited to this spatial scale as the lensing analysis 
was performed on HST-WFPC2 data with pointings at cluster centers. The remaining 80- 
90% of the cluster mass is consistent with being smoothly distributed (in lumps with mass 
M < 10 10 M o ). 

In Fig. 1, we show the mass function retrieved from galaxy-galaxy lensing for each of the 
five HST clusters. There is a low-mass cut-off in the observed clump spectrum, at around 
1O 11 M , which is due to observational limitations. The mass resolution of this technique 
is limited by the depth and field of view of the Wide Field Planetary Camera (WFPC2) 
aboard the Hubble Space Telescope, by the number of background galaxies per foreground 
lens, and the reliability with which shapes can be measured for the faintest background 
galaxies in the HST image of the central region. Unfortunately, this limits the number of 
reliably determined lumps per rich cluster to about 40, implying that only the massive end of 
the clump spectrum corresponding to the brightest cluster galaxies can be probed. However, 
despite the low number statistics, the individual cluster measurements show a marked rise 
in substructure abundance towards lower mass scales. This becomes particularly apparent 
once the clusters are stacked, as we have done in the lower right panel of Fig. 1. 



3. Cluster substructure in high resolution ACDM simulations 

We analyze substructure in high-resolution dark matter simulations of clusters formed 
in the ACDM model. Our set of simulations is taken from the study carried out by Springel 
et al. (2001a) of a single rich cluster of mass 8.4 x lO 14 /i _1 M , simulated in 4 steps of 
ever increasing resolution using the parallel tree-code GADGET (Springel et al. 2001b). In 
these simulations (referred to as 'SI' to 'S4'), the particle mass resolution increases from 
6.87 x 10 9 to 4.68 x 1O 7 /i~ 1 M , corresponding to 1.3 x 10 5 particles up to about 20 million 
within the virial radius, making the highest resolution simulation in this series one of the 
best resolved simulations of a single rich cluster carried out to date. This simulated cluster 
has a comparable central density and mass within the inner 500 h^ 1 kpc as the massive HST 
cluster-lenses studied here. 

Self-bound gravitational substructure is found with the algorithm SUBFIND (Springel 
et al. 2001a). It starts by determining locally over-dense dark matter substructure candidates 
in a fully adaptive fashion, and then subjects each of them to a gravitational unbinding 
procedure, such that a catalog of self-bound dark matter substructures results. In the lowest 
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resolution simulation SI of the series, 118 substructures can be detected, a number that 
increases to more than 4600 in the high-resolution simulation S4. Note however that these 
additionally resolved small substructures are all of ever lower mass; already the low resolution 
simulation captures the correct number of massive satellites, with (at least on average) the 
correct mass. 

The measured differential mass function of substructures appears to be a power law, 
dN/dm oc m~ a , with slope close to a ~ —1.8 (Springel et al. 2001a; Helmi et al. 2002; 
De Lucia et al. 2004). Interestingly, this is very close to the result by Lee (2004) who 
has attempted an analytic calculation of the sub-halo mass function based on modeling the 
complex dynamical history of galaxies in the cluster with a parameterized model to account 
for the effects of global tidal truncation (see also Taylor & Babul 2001). Given that we 
have also demonstrated that the spatial extents of substructures inferred from galaxy-galaxy 
lensing are consistent with the tidal stripping hypothesis, this lends support to the validity 
of a simple tidal-limit approximation. Similar to the lensing result, both the analytic model 
of Lee (2004) and the numerical simulations find that only about 10% of the sub-halo mass 
is bound in substructures, most of it in a handful of most massive sub-halos which dominate 
the cumulative mass in the substructures. 

We note that numerical studies by De Lucia et al. (2004) find that the substructure 
mass function depends only weakly on the properties of the parent halo mass. Also, the 
mass fraction in substructure is relatively insensitive to the tilt and overall normalization of 
the primordial power spectrum (Zentner & Bullock 2003). Only for radically altered CDM 
models, for example by truncating small-scale power, Zentner & Bullock find that their 
models yield projected substructure mass fractions that are lower than the estimates from 
strong lensing. 

4. Comparison 

We now compare the galaxy-galaxy lensing results with the substructure mass function 
obtained from the N-body simulations. In Fig. 2, we show histograms for the distribution 
of substructure masses in the four simulations of Springel et al. (2001a) and contrast them 
individually with the stacked result for the clusters A2218, A2390, and C10054. Since the 
number of substructures in the observable mass range is quite small, we expect large system- 
to-system variations between different clusters. The stacking of clusters of similar mass that 
we applied here reduces the associated scatter somewhat. 

Comparing the lensing result with the four simulated clusters (which also show some 
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numerical scatter amongst each other), we find broad general agreement in the mass range 
10 11 — 10 12 5 M Q , both in the amplitude and the slope of the substructure mass function. 
Given that no free parameter or scaling has been applied to obtain this match, the agreement 
is in fact remarkable. 

However unlike the observations, the simulations show no low-mass cut-off. If this cut- 
off is entirely due to resolution effects, which appears plausible, the interpretation is that the 
observations only see the 'tip of the iceberg' of the substructure distribution as far as their 
number is concerned, even though they detect most of the mass in substructures. Fig. 2 
gives also a hint that the simulations may systematically over-predict the mass of the most 
massive satellites, but we caution that such an effect could also be caused a systematic effect 
in the mass estimate based on the lensing technique, for example if the massive substructures 
preferentially correspond to halos that have fallen in most recently. 

5. Discussion 

On comparing the clump mass spectrum obtained for galaxy clusters from high-resolution 
N-body simulations of ACDM models to those obtained from galaxy-galaxy lensing in HST 
cluster lenses, we find excellent agreement. Despite the fact that the lensing analysis assumes 
that mass traces light and is only sensitive to a restricted mass range, it is clear that there 
is no substructure problem in CDM on mass scales spanning 10 11 — 1O 125 M . This is in 
sharp contrast to the situation on galactic scales, where the paucity of observed satellites 
when compared with the abundant dark matter substructure predicted by simulations has 
been characterized as a crises for CDM. While the severity of this problem has probably 
been overstated initially - in fact it may have partially gone away by now (Stoehr et al. 
2002) - it is clear that CDM predicts a rich spectrum of dark matter substructure extending 
to very small masses, both for clusters and galactic halos. In clusters, up to several hun- 
dred of the most massive substructures can be directly identified with the luminous cluster 
galaxies, and semi-analytic models of galaxy formation show that this association leads to 
highly successful models for the population of cluster galaxies (Springel et al. 2001a). On 
the other hand, on galactic scales, a much larger number of truly dark satellites must be 
present according to the N-body models. Gravitational lensing is probably our best bet to 
detect such structures of small mass, and in fact, on these small scales, the observed flux 
anomalies in multiply-imaged quasar systems (Mao & Schneider 1998; Chiba 2002; Dalai & 
Kochanek 2002; Metcalf & Madau 2001; Metcalf & Zhao 2002; Mao et al. 2004) have been 
interpreted as evidence for significant substructure in the mass range 10 4 < M < 10 8 M Q . 

In this letter, we have shown that the clump mass function independently determined 



- 8- 



from lensing (a technique that is unaffected by the dynamical state of the cluster) is in 
excellent agreement with that obtained in high resolution cosmological N-body simulations 
of clusters of galaxies in the ACDM model. 

PN thanks her collaborators on the HST cluster-lenses project: Jean-Paul Kneib, Ian 
Smail and Richard Ellis. 
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Fig. 1. — The substructure mass function retrieved from galaxy-galaxy lensing for 5 massive 
HST cluster-lenses. This mass spectrum is constituted by the dark halos associated with the 
40 brightest galaxies in these clusters. The panel on the bottom right shows the coadded 
and averaged clump spectrum of clusters A2218, A2390, and C10054. 




Fig. 2. — Comparison of the substructure mass function of the stacked clusters A2210, 
A2390, and C10054 (shaded histogram) with N-body simulations of a rich cluster of galaxies 
(filled histogram), carried out at four different numerical resolutions. 



